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ABSTRACT 


A new second-order theory was developed for predicting the service lives of aerospace structural 
components. The predictions based on this new theory were compared with those based on the Ko 
first-order theory and the classical theory of service life predictions. The new theory gives very accu- 
rate service life predictions. An equivalent constant-amplitude stress cycles method was proposed for 
representing the random load spectrum for crack growth calculations. This method predicts the most 
conservative service life. The proposed use of minimum detectable crack size, instead of proof load 
established crack size as an initial crack size for crack growth calculations, could give more a realistic 
service life. 


INTRODUCTION 

Service life (or the number of flights permitted within the proof load interval) of an aerospace ve- 
hicle is governed by the individual service lives of critical structural components (for example, NASA 
B-52 carrier aircraft air-launching-system hooks (ref. 1)). The following procedure describes the con- 
ventional method of estimating the service life of critical structural components. 

First, a proof load test (which covers all critical flight maneuver loading conditions) is conducted on 
the critical structural components to load those components up to design-limit load levels. If the proof 
load test should cause destruction on certain components, then those failed components are immediately 
replaced, and another proof load test is repeated. If all the structural components should survive the 
new proof load test, then fracture mechanics is used to theoretically determine the initial fictitious 
(nonexistent) crack size, Oc. (critical or incipient crack size associated with the proof load level) at 
the critical stress point (the location of which is determined from stress analysis) of each structural 
component. Then, the service life for each structural component is estimated from the amount of crack 
growth permitted for the initial “fictitious” crack a % to grow and reach the operational limit crack size 
a°(a° > a?) which is also fictitious, and is calculated from fracture mechanics based on the operational 
peak load, which is much lower than the proof load. Thus, the amount of available crack growth a® — o? c 
will determine the number of flights available for operation until the next proof load test. 

Very often the initial “fictitious” crack size a? established through the proof load test is much 
larger than the minimum observable crack size (a 0 ). Therefore, the service life predicted using the 
proof load established initial crack size aP c could be too conservative and unrealistically short compared 
with the service life predicted using the minimum observable crack size ( a Q ) as the initial crack size in 
the service life calculations. 

Sometimes, it might be more convenient to relate the amount of crack growth caused by random 
load spectrum to the number of equivalent constant-amplitude stress cycles, and then estimate the service 
life by simply counting the available constant-amplitude stress cycles. 

A second-order theory for calculating the service life (number of remaining flights) for a given 
available crack growth (refs. 1—3) is presented in this report. In addition, a discussion of the previ- 
ously mentioned equivalent constant-amplitude stress cycling method to estimate the service lives of 
aerostructural components is also presented. 
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crack location parameter 

depth of semielliptic surface crack, in. 

operational limit crack size, in. 

initial fictitious crack size established by proof load test, in. 
minimum detectable crack size, in. 
crack size at the end of the first flight, in. 

material constant in Walker crack growth-rate equation, c ^ e ^si\/im)~ m 

half length of surface crack, in. 

complete elliptic function of the second kind 

number of remaining flights calculated based on the first-flight data 

number of remaining flights calculated from the equivalent constant-amplitude stress 
cycles method 

number of remaining flights calculated from the first-order theory 

number of remaining flights calculated from the second-order theory 

fraction of limit load (proof load) 

mode / critical stress intensity factor, ksi\/Tm 

mode I stress intensity factor associated with cr max 

modulus of elliptic function 

number of flights 

flaw magnification factor 

Walker exponent associated with stress amplitude 
number of constant-amplitude stress cycles 

maximum number of constant-amplitude stress cycles for service life 

number of random stress cycles 

Walker exponent associated with stress ratio 

surface flaw and plasticity factor 

stress ratio, R = &mo.x 

solid rocket booster drop test vehicle 

front hook vertical load, lb 

left rear hook vertical load, lb 

right rear hook vertical load, lb 

amount of crack growth induced by the first flight, in. 

number of equivalent constant-amplitude stress cycles consumed during the £th flight 
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Sai crack growth increment induced by the ith half cycle, in. 

a* proof load established stress at the critical stress point, ksi 

a j tensile stress at stress point j 

a max maximum stress of a stress cycle, ksi 

a min minimum stress of a stress cycle, ksi 

<7 S mean stress (or static stress) of a stress cycle, ksi 

an ultimate stress, ksi 

ay yield stress, ksi 

r ultimate shear stress, ksi 

4> angular coordinate for semielliptic surface crack, rad 

( )j quantity associated with ith half stress cycle 

SERVICE LIFE 


Conventional Method 

If Aoj is the amount of crack growth induced by the first flight, then the conventional method 
predicts the number of remaining flights F\ (service life) based on the following equation (refs. 1-3). 


El = 


a r 


P 

ac 


Aai 


( 1 ) 


where aP c and a° c are calculated respectively from (refs. 1-3) 


oP c = ~ 

" 7 r 


K 


*c 


AM k o * 


Q 


a r = — 

C 7T 


Kir 


1 2 


(2) 


(3) 


AM K fa 

where a* and /a* (/ < 1) are respectively, the proof load induced stress (limit stress) and the op- 
erational peak stress at the critical stress point; A is the crack location parameter ( A = 1.00 for the 
through thickness crack, A = 1.12 for both the surface and the edge cracks)(refs. 1-3); M K is the flaw 
magnification factor {Mk = 1.0 for very shallow surface cracks, Mk =1.6 when the depth of the crack 
approaches the thickness of the plate); Kj c is the critical stress intensity factor, and Q is the surface 
flow shape and plasticity factor of a surface crack. If the surface crack is semielliptic in shape, then Q 
is expressed as (see fig. 1 and refs. 1-3): 


Q = [E(k)} 2 -02n(^J 


(4) 


where ay is the yield stress, and E(k) is the complete elliptic function of the second kind defined as 


E{k ) = v/l - fc 2 sin 2 </> d(j) 
Jo 


(5) 
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where <f> is the angular coordinate for a semielliptic surface crack (see fig. 1) and the modulus k of the 
elliptic function is defined as 

k = 

where a and c are respectively the depth and the half length of a semielliptic surface crack. 

Before flight, the actual amount of crack growth Aai (eq. (1)) for the first flight is unknown. The 
way to estimate Aai, before the actual flight, is to perform transient dynamic analysis of the flight 
vehicle under specified severe maneuvers such as landing, braking, ground turns, flight in severe buffet 
and turbulence, etc. Actual ground maneuvering of the aircraft can be conducted (for example, taxi 
runs on straight or curved paths) and generate an actual loading spectrum for each critical component 
for a short period. Then, the loading spectrum is expanded (extrapolated) to cover the duration of one 
flight. For large flexible aircraft, such as the B-52 carrier airplane, the ground maneuver could produce 
a more severe loading spectrum than that of the actual steady flight. If equation (1) predicts a sufficient 
number of flights available based on Aai, calculated from the ground maneuver, one may feel confident 
to actually conduct one flight to obtain the actual value of Aai. 

Calculation of Crack Growth 

The crack growth caused by the random stress cycling of the first flight may be calculated by using 
the half-cycle theory (refs. 1-3). The half-cycle theory states that the damage, or crack growth, caused 
by each half cycle (either increasing or decreasing load) of the random load spectrum is assumed to 
equal one half of the damage caused by a full cycle of the constant-amplitude load spectrum of the same 
loading magnitude. Thus, the total damage done by the random load spectrum will be the summation of 
the microdamages caused by the individual half waves of different loading magnitudes. Figure 3 gives 
graphical illustrations of the half-cycle theory (refs. 1-3). 

Thus, the crack growth Aaj caused by the first flight may be calculated as 

2N n 

Aaj = a! - aP c = £ tfa* (7) 

1 = 1 

where a\ is the crack size at the end of the first flight; Nr is the total number of random cycles induced 
by the first flight; and <5aj is the crack growth increment induced by the ith half cycle, calculated from 
(refs. 1-3) 

Sai^^KKmaxhni-Ri) 71 ( 8 ) 

which is obtained from the Walker equation (refs. 1-3) 

A = C(K max r( 1 - R) n (9) 

by setting da = Sa^dN = j,K m ax — ( K ma x)i and R = R^ for the zth half cycle of random load 
spectrum. In equations (8) and (9), C, m, and n are the material constants and ( K ma x)i and R{ are, 
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respectively, the maximum stress intensity factor and the stress ratio associated with the zth half cycle 
and are given by 

( Kmax)i — A Mk (cr max) i 

r _ = (grmn li (11) 

\&max)i 

where (cr max )j and {cr m i n )i are respectively, the maximum and the minimum stresses of the zth half 
cycle (see random stress cycles in fig. 2(a)); Qj_i is the crack size at the end of the (i— l)th half cycle. 


“Minimum-Crack” Method 


The proof load established initial crack sizes (q?) are used only for establishing a “baseline” for the 
aircraft structural component already in service. If the critical stress-point areas can be easily inspected, 
then the minimum detectable crack size a 0 could be used as an initial crack size. Certainly aircraft with 
zero flight hours can use a 0 as an initial crack size with no reliance on proof load requirement. 


If the minimum detectable surface crack size a Q (crack depth) turned out to be much smaller than 
the proof load established intial crack size q§ of certain critical structural components (for example, the 
B-52 hooks which are inspectable) then one can use a Q instead of a? as an initial “fictitious” crack size 
for the crack growth calculations. Thus, equation (1) may be modified as 


_ qg-a 0 
Fl A a x 


( 12 ) 


where Aqj is much smaller than Aqj appearing in equation (1) because of smaller initial crack size. 


If equation (12) is used in the calculations of service life, the initial crack size a 0 for all subse- 
quent flights will remain the same provided interflight crack detection inspection is conducted. Clearly 
equation (12) will give much longer service life than does equation (1). 


SECOND-ORDER THEORY 

The conventional equation (eq. (1)) for service life prediction is based on the assumption that the 
amount of crack growth during each flight for all subsequent flights remains the same as the crack growth 
Aqj, caused by the first flight. In reality, the amount of crack growth during each flight will increase 
steadily with the number of flights accumulated because the initial crack size for the subsequent flight 
will increase gradually. The new equation for service life prediction will account for the nonuniform 
crack growth effect. If Aq^ is the amount of crack growth induced by the random load spectrum of 
the Ah (i = 1, 2, 3, . . .) flight, and if A N t is the number of cycles of an equivalent constant-amplitude 
load spectrum which also induce a crack growth equal to A q^, then the Walker equation (eq. (9)) may 
be used to relate A to A a,j> as 


Aq^ — C 


^AMxOmax 



m 


(1 


R) n (ae-i)^AN e 


(13) 
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in which the following expression was used 


Kmax — AMKamax^ — Q~^ 


where a^_j is the crack size at the end of the (£— l)th flight. 


(14) 


For simplicity, if we assume that the equivalent constant-amplitude load spectra for all flights are 
identical (that is, cr-max » -R> and A/V^ remain the same), then equation (13) could be used to establish 

the following crack growth ratios and expand them in terms of ^ L up to second-order terms assuming 

that — tjI- is small (i.e., < 1): 

a c a c 
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If the available crack growth a° c — afc can allow F\ number of flights, then one can write: 

F\ terms 


a o a P A a, i + Ad2 ""b Aa 3 *b • • • “b A a,£ + . . . Aa 


F i 


Aai 


Aoi 


( 20 ) 


where the left-hand side of this equation is F\ according to equation (1). 
Substitution of equations (15)— (19) into equation (20) yields: 

F i terms 

F\ = r+TTT+T+TTTTT 


+ 


+ 


(Fj-l) terms 


' m A a 


1 


0 + 1 + 2 + 3 + . . . + (t - 1) + • ■ • + (Fi - 1) J 

(F\ — 1) terms 


2 dr 


l { m A d\ 

1_ 2 lY^T 


0 2 + l 2 + 2 2 + 3 2 + . . . + (i - l) 2 + • • • + (Fi - l) 2 
which becomes, after summations of numbers: 


(’ 4 ) 


1 \ ( m Aai 


+ 


1 - - . ( m Aai 

F^F l + -F l( F-l) 


^ 1 ( m Aaj 


2\2 a p c 

1 - - / - 1 \ / 1 \ fm Aa] 

+3 Fm - 1) (F, - 5 ) (*--) 


+ . .. 


( 21 ) 


( 22 ) 


After grouping terms and rearrangement, equation (22) may be written in the following standard form 
of cubic equation: 

F\ > + pF \ 2 + qF\ + r » 0 


(23) 


where 


3 (, 3 

^2(^T)( 1 "2 m + 2 A^. 


q = \v + ^ T 


m _ 2 Q c + / a c \ 2 

2 Aai m \^ a l J 




\2Fi 


p \ 2 

Oc \ 


m(m — 1) \Aai 
The real root of equation (23) is then given by 


( 24 ) 


( 25 ) 


( 26 ) 


F\ = B + D - 


( 27 ) 
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where 


where 


i>4 

(28) 

Q = ^(3 q-p 2 ) 

(29) 

P = ^( 2 P 3 “ + 27r) 

(30) 


In equation (21), if the second-order terms are neglected, then the number of flights F x predicted based 
on the first-order expansion will be 


r; 2a.c 

F ' = ^[ 



m ^F, - l) 

O-C J 


which has already been published in references 1 and 2. 


(31) 


EQUIVALENT CONSTANT-AMPLITUDE STRESS CYCLES 

In this section, we attempt to consider the crack growth A a lf (see eq. (7)) caused by random 
stress cycling as if it were caused by equivalent constant-amplitude stress cycling (see fig. 2(b)). If the 
constant-amplitude load spectrum is cycling about the mean (static) stress a s , then the maximum stress 
a max and the minimum stress a m j n of the equivalent constant-amplitude stress cycle is related through 


° min — 2(7 s &max 

^ a max is a fraction of the limit stress a* (proof load induced stress), namely. 


(32) 


then one can write 


&max — /<7 ; 


/< 1 


K, 


max 


= AM K fo* 



(33) 

(34) 


R = 


a min 

<7 max 


2ct 5 

7^ 


(35) 


If the initial crack size is a§, then the number of the equivalent constant-amplitude stress cycles, N, 
required to extend the initial crack size a c up to an arbitrary size “a” could be obtained by integrating 
equation (9) as 


N = 


(* 


a 


1-f - 


-K ) 1 


m 

T 


m 

T 


) c 


AM K fa* 



m 




(36) 
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From equation (36), the number of equivalent constant-amplitude stress cycles N c available for opera- 
tion, and the number of equivalent constant-amplitude stress cycles consumed by the first flight, N lt to 
cause the same damage A a\ as does the random stress cycling, may be obtained by setting a — a° and 
a = aj respectively. 


If every flight after the first flight consumes the same number of equivalent constant-amplitude 
stress cycles as does the first flight, then the service life F f may be calculated from 




N c 
N i 


1 -^ 


- (a£) 


where equation (36) is used, and a\ — oP c + Aa^ is to be calculated from equation (7). 


(37) 


When the random load spectrum is converted into an equivalent constant-amplitude load spectrum, 
for service life estimates one can merely count the number of the equivalent constant-amplitude stress 
cycles. If N\ is consumed for the first flight, then N\ is used to calculate the remaining flights (eq. (37)) 
with accuracy, because the amount of the equivalent constant-amplitude stress cycles consumed for each 
of the subsequent flights is theoretically identical. 


NUMERICAL EXAMPLE 

The numerical example chosen for the calculations of service life of aerospace structural components 
is the NASA B-52 carrier aircraft air-launching system (pylon) hooks as shown in figures 4 and 5. 
Figure 5 also shows the locations of the critical stress points and the stress— hook-load relationships 
determined from stress analysis (ref. 4). Through one front hook and two rear hooks of the pylon of 
the B-52 aircraft, a winged Pegasus® rocket (approximately 44,629 lb) will be carried to high altitude 
(approximately 40,000 ft) to release it for firing into orbit. The B-52 aircraft was used earlier to carry 
heavy stores such as the X-15 air-launched rocket plane (57,250 lb) and the space shuttle solid rocket 
booster drop test vehicle (SRB-DTV, 49,000 lb). The data accumulated for those vehicles may be used 
to estimate the preflight service lives of the three hooks when the store is the Pegasus rocket, because 
of weight proximity. 


Input Numerical Values 

^ Assuming that the surface cracks (initial and after growth) are semicircular in shape (that is, 
2c = l)' anc * ^ at the stress level at the critical stress point of each hook reached the yield point (the 
hooks are designed to carry yielding zones), then from figure 1 one obtains Q = 2.265 for ^ In 
the crack growth calculations A = 1.12 and M x =1.0 were used. Other numerical values used in the 
crack growth calculations are given in table 1. 


® Pegasus is a registered trademark of Orbital Sciences Corp., Fairfax, Virginia. 
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Table 1. Material properties for B-52 pylon hooks. 


Part 



oy. 

T, 

Ki c > 

c, 



name 

Material 

ksi 

ksi 

ksi 

ksi Vxm 

sdfe (ksi v/taT) m 

TO 

n 

Front hook 

Inconel 718® alloy 

175 

145 

135 

125 

9.220 xlO" 12 

3.60 

2.16 

Left rear hook 

AMAX MP35N* alloy 

250 

235 

141 

124 

2.944 x 10 -11 

3.24 

1.69 

Right rear hook 

AMAX MP35N* alloy 

250 

235 

141 

124 

2.944 x 10 -11 

3.24 

1.69 


® Inconel 7 18 is a registered trademark of Huntington Alloy Products Division, International Nickel Company, Huntington, 
West Virginia. 

*AMAX MP35N is a trademark of SPS Technologies, Inc., Jenkintown, Pennsylvania. 


Using the numerical values given in table 1, and the proof hook loads Va = 36,500 lb, V BL = 
Vbr — 57,819 lb, the initial crack sizes a? c and the operational limit crack sizes a° (for / = 0.6) may 
be calculated respectively from equations (2) and (3) as: 


Table 2. Proof and operational limit crack sizes. 


Part 

name 

Proof hook 
load, lb 

p ■ 

etc. m. 

a° c , in .(/ = 0.6)* 

^ (/ = 0.6)* 
dr 

Front hook, Va 

36,500 

0.1247 

0.3465 

0.01814 

Left rear hook, V BL 

57,819 

0.0774 

0.2151 

0.00761 

Right rear hook, V BB 

57,819 

0.0774 

0.2151 

0.00761 


*/ = 0.6 was the average operational peak load level for the case of SRB-DTV. 


The present-day crack detection techniques could detect a surface crack approximately 0.02 in. 
long. Thus, if the surface crack is semicircular in shape, then the minimum detectable crack depth is 
approximately 0.01 in. Clearly, this value is far less than the proof load established initial crack sizes 
listed in table 2. In the service life calculation using the “Minimum-Crack” method, the initial crack 
size a 0 will be taken as a 0 = 0.01 in. 


Crack Growth Ratios 


Figure 6 shows crack growth ratio plotted as a function of the number of flights £ for the 
front hook using the first- and the second-order theories. The first-order theory (eq. (31)) gives linear 
increase of with increasing £. However, the second-order theory (eq. (23)) gives nonlinear curve 

with lower values of in the region 1 < l < 41, beyond that the second-order theory predicts much 

higher values of than the first-order theory. The conventional theory gives a simple horizontal line 


Aa 
A a 


* = 1 . 
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Figure 7 shows similar plots for the rear hook. Unlike the front hook, the second-order curve does 


not intersect with the first-order curve until £ ~ 120 because of different values of m and 


Aai 


Remaining Flights 

Table 3 lists the numbers of remaining flights F\ and F\, for the front and rear hooks, calculated 
respectively from the first- and the second-order theories, compared with the corresponding number of 
flights F\ based on the conventional method (assuming equal amount of crack growth for all subsequent 
flights (ref. 1)). It is clear that the conventional theory exceedingly overpredicts the predicted service 
life. Values in table 3 are plotted in figures 8 and 9 respectively for the front and rear hooks for easy 
visualization of the curves for F \ , Fj, and F\ . 


Table 3. Number of remaining flights calculated from first- and second-order theories compared to the 
conventional method. 


F i 

Front hook 

Rear hook 

h 

F i 

Fi 

F , 

1 

i 

1 

1 

1 

50 

33 

30 

40 

39 

100 

53 

45 

70 

65 

150 

70 

56 

95 

85 

200 

84 

65 

116 

101 

250 

97 

72 

136 

116 

300 

108 

79 

154 

128 

350 

119 

85 

171 

139 

400 

129 

90 

186 

149 

450 

138 

95 

201 

159 

500 

147 

99 

215 

167 

550 

155 

103 

228 

175 

600 

164 

107 

241 

183 

650 

171 

111 

254 

190 

700 

179 

114 

265 

197 

750 

186 

118 

277 

204 

800 

193 

121 

288 

210 

850 

200 

124 

299 

216 

900 

206 

127 

310 

222 

950 

213 

130 

320 

227 

1000 

219 

132 

330 

233 


For the given values of m and — p *- for the B-52 front and rear hooks (tables 1 and 2), the predicted 

a c 

service life for the two hooks based on the conventional theory and the first- and the second-order theories 
are presented in table 4: 
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Table 4. B-52 hooks service lives based on conventional, first- and second-order theories. 



F i 

A 

A 

Front hook 

98 

52 

45 

Rear hook 

234 

130 

110 


Clearly the conventional theory gives too optimistic a service life prediction. The service life of 
the front hook is somewhat shorter than that of the rear hook. 

Figures 10 and 11 show the crack growth curves for the front and rear hooks, respectively, based 
on the equivalent constant-amplitude stress cycle method. In each figure, the upper curve is the plot of 
equation (37) (that is, initial crack size is the proof load established crack size a?). The lower curve is 
the plot of equation (37) with etc replaced with a 0 as an initial crack size. Using a 0 as an initial crack 
size instead of etc, the service lives of the hooks are greatly enhanced. Table 5 summarizes the results 
of the equivalent constant-amplitude stress cycle method. Predictions from the conventional theory are 
also shown for comparison. 

Table 5. B-52 hooks service lives based on the equivalent constant-amplitude stress cycles method, 
(/= 0 . 6 ). 




Initial crack size 




P 

a c do 

P 

a c do 

V 

a c do 

a p c 


Stress cycles 

Remaining flights 



Ni 

N c 

A* 

A 

Front hook 

67.08 67.08 

2623 3329 

39 496 

98 

Rear hook 

14.88 14.88 

1490 9602 

100 645 

234 


Notice that the number of equivalent constant-amplitude stress cycles for cycling at 60 percent 
(/ = 0-6) of limit stress consumed during each flight (N\) are relatively low and are independent of 
the initial crack size. For the initial crack size of a p , the equivalent constant-amplitude stress cycles 
method predicts the most conservative service life as compared with other theories (tables 4 and 5). 

Accuracies of Expansions 

To check the accuracies of the first- and the second-order theories, values of the crack growth rate 
( ec l- (16)) were calculated from the exact expression for (before expansion of eq. (16)), and 
from the first- and the second-order expansions (eq. (16)). The results are tabulated in table 6. 
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Table 6. Accuracy of expansion for ^2. 

A 


A°2 

Aa i 

Front hook 

Rear hook 

Exact expression 
(>^) f 

1.032887789 

1.012357255 

First-order expansion 

, * 

He 

1 + ^ Aa i 

: 1. 032 j65 1163 

■ 1.0123 -28200 

a c 

Second-order expansion 

* 

% 

! 1.0123572 i 82 

(eq. (16)) 

il. 03288 j 8074 


j accurate digits 


This table shows that the second-order expansion gives very accurate values for , thus, the 
service life predicted from the second-order theory is quite reliable. 

CONCLUDING REMARKS 

The second-order theory for predicting the service life of aerospace structural components was 
presented. The service life predicted from the second-order theory was compared with those predicted 
from the previously developed first-order theory and the conventional method (constant amount of crack 
growth for all subsequent flights) of service life predictions. The second-order theory (based on the 
second-order expansion of the crack growth rate) could give reasonably accurate values of crack growth 
rate compared with the exact values. 

The new equivalent constant-amplitude stress cycles method was proposed. This method gave the 
most conservative service life predictions. The use of minimum detectable crack size, instead of proof 
load established crack size as an initial crack size, could give a more realistically longer service life. 

Dryden Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, August 18, 1991 
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Figure 1. Surface flaw shape and plasticity factor for semielliptic surface crack. 
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(a) Random stress cycles. 



(b) Constant-amplitude stress cycles 


Figure 2. Random stress cycles and equivalent constant- 
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implitude stress cycles. 


'-■jjj = C(AK) m (1 -R) n_ 





B-52 wing 



Figure 4. B-52 aircraft pylon carrying winged Pegasus rocket. 







Figure 6. Crack growth ratio as a function of number of flights based on the first- and second-order 
theories; front hook. 
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Figure 7. Crack growth ratio as a function of number of flights based on the first- and second-order 
theories; rear hook. 


20 





Figure 8. Comparison of number of remaining flights predicted from conventional theory and the first- 
and second-order theories; front hook. 
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Figure 9. Comparison of number of remaining flights predicted from conventional theory and the first- 
and second-order theories; rear hook. 
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Figure 10. Comparison of service lives of front hook based on different initial crack sizes. 



Figure 11. Comparison of service lives of rear hook based on different initial crack sizes. 
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